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We examine the revival features in wave packet dynamics of a particle confined in a finite square 
well potential. The possibility of tunneling modifies the revival pattern as compared to an infinite 
square well potential. We study the dependence of the revival times on the depth of the square 
well and predict the existence of superrevivals. The nature of these superrevivals is compared with 
similar features seen in the dynamics of wavepackets in an anharmonic oscillator potential. 
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I. INTRODUCTION 

The wavepacket in quantum mechanics is often viewed 
as the most 'nearly classical' state and is known to exhibit 
many striking classical properties. However, its inherent 
quantum nature also causes it to exhibit many quantum 
mechanical features. It is known that in certain 
nonlinear quantum systems the dynamics of the wave 
packet incorporates quantum interference effects which 
cause it to undergo a sequence of collapses and revivals 
and in the course of its evolution the wavepacket period- 
ically breaks up and reconstitutes its original form. At 
intermediate times the wavepacket gathers into a series 
of subsidiary wavepackets called fractional revivals ||. 
Quantum revivals are often a manifestation of the fact 
that the time evolution of the wavepacket is driven by a 
discrete eigenvalue spectrum and revival features depend 
on the way the eigenenergies of the quantum system de- 
pend on the quantum number, n |||]]. Recently there 
has been a lot of interest in the theoretical and experi- 
mental study of quantum revivals in a variety of nonlin- 
ear systems like that of Rydberg atom wavepackets [||, 
molecular wavepackets ||, wave packets in semiconduc- 
tor quantum wells Q etc. Most experiments dealing with 
quantum revivals till date focus on atomic 0] and molec- 
ular systems ||, photon cavity systems and ions in 
traps ||. The simplest class of systems for which one 
would see fractional and full revivals are those for which 
the energy spectrum goes as n 2 , e.g., the infinite square 
well potential and the rigid rotator [Q . For systems with 
a nonquadratic dependence on n, one can see a new se- 
quence of collapse and revivals which are different from 
the usual fractional revivals. Often these sequences cul- 
minate with the wavepacket resemblingits original form 
more closely. These are superrevivals [QJ|. The revival 
patterns of Rydberg wavepackets, where the energy spec- 
trum is nonquadratic in n, have been seen to exhibit such 
superrevivals jjj. 

An ideal model system to illustrate the fractional re- 
vival phenomenon is the infinite square well potential 
since simple analytical forms of the eigenvalues and the 



eigenf unctions allow for an easy analysis of the time evo- 
lution of any initial state. This system has been recently 
studied in great detail in the context of fractional re- 
vivals Jl(j,[ll| as well as for the space-time structures O] 
that appear in the dynamical evolution. An experimen- 
tal realization of the predictions of the fractional revival 
phenomenon in the infinite well is most likely to be in 
semiconductor quantum well systems. In reality, how- 
ever, it is rather impossible to find a physical system 
that creates a truly infinite confining potential. It be- 
comes crucial, therefore, to to study the problem of re- 
vivals in more realistic, physical systems which are bet- 
ter described by finite well potentials rather than infi- 
nite potentials. The possibility of the occurrence of both 
revivals and superrevivals in a finite potential well has 
been pointed out before pfj| though no explicit study 
has been done till now. From a more fundamental point 
of view, the motivation to study the finite well system is 
to get a greater insight into understanding the 'classical 
limit'. During the course of its evolution an initially lo- 
calized wavepacket appears at certain times as a linear 
superposition of spatially separated copies of itself, i.e., 
in Schrodinger cat-like states. The tremendous progress 
in experiments involving semiconductor systems, e.g, the 
observation of quantum beats in quantum wells 31 and 
Bloch oscillations in semiconductor superlattices | K| in- 
creases the prospects of generating, detecting and study- 
ing such cat-like states among other features like revivals, 
thus providing further motivations to get a better theo- 
retical understanding of these systems. 

In this paper we focus on the revival features of 
wavepackets in finite square well potentials and com- 
pare it with the corresponding case for the infinite well. 
For our purpose we have only concerned ourselves with 
the bound eigenstates of the finite square well potential. 
According to quantum mechanics, the wavefunction has 
nonvanishing values in the classically forbidden regions 
thus giving a nonvanishing probability for the particle 
being outside the well. One thus expects a difference 
in the wavepacket dynamics as compared to the infinite 
well case. We show that the existence of the possibil- 
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ity of tunneling modifies the revival pattern as compared 
to the infinite square well potential. In particular it al- 
lows for revivals and superrevivals. We show that usual 
revival times are now longer compared to that of the in- 
finite well and depend on the well depth. The paper is 
organized as follows. In Section II we begin by introduc- 
ing the concept of fractional revivals and review briefly 
the finite well problem. We then present our numerical 
results and confirm the approximate formula of Barker et 
al. [ fl4[ from which we deduce an approximate analytical 
expression for revival times in finite square well poten- 
tials. In Section III we present another key result of this 
paper which is the existence of superrevivals during the 
time evolution of an initial Gaussian wavepacket. The 
absence of an exact analytical expression for the energy 
spectrum for the finite potential well makes it difficult to 
give an estimate of the superrevival times. However, we 
get some insight into the nature of superrevivals in the 
well by comparing it with superrevivals in the dynamics 
of Gaussian wavepackets in an anharmonic oscillator po- 
tential for which we have the analytical results. Finally, 
we summarize our results in Section IV. 



II. DEPENDENCE OF REVIVAL TIMES ON 
WELL-STRENGTH 

Consider the time evolution of a particle initially in 
state tp(x,0) in a potential: 



tp(x,t) = y^c„0»(x)e" 



■iE n t/h 



(1) 



etc. which are generally termed in the literature as the 
classical, revival, superrevival times and so on [Q. One 
can easily see that re- writing the time evolution (1) in 
terms of these time scales shows how they govern the 
time evolution of ?/;(a;,0). The time scales in turn are 
controlled by the dependence of the energy on the quan- 
tum number n. For the simple case of the infinite well, 
the quantized energy levels are exactly quadratic in n: 



7T 2 ftV 

2mL 2 ' 



(4) 



where m is the mass of the particle and L is the length of 
the well. Corresponding to this, one has the time scales 
T c i = 2mL 2 /irfin , T rv — AmL 2 /itfi while T sr = oo. The 
time evolution (1) can be re-written as: 



ip(x,t) = } j c n (f} n {x)e 



-2Tri(t/T r v)n 2 



(5) 



The expansion (5) includes both odd and even parity 
states. It is easy to see that the wavefunction regains its 
original form i.e., shows full revivals, whenever t equals 
some multiple of T rv , and shows fractional revivals when- 
ever t is equal to some rational fraction of T rv , e.g, T r „/4, 
but no higher order effects like superrevivals are seen in 
the dynamics of the wavepacket for the infinite square 
well potential JT~0| ] . 

We now turn to the problem of the finite square well 
potential which is the focus of this paper. The one- 
dimensional finite square well potential for a well of 
length L can be described as: 



where n is the quantum number and 4> n (x) and E n are the 
energy eigenstates and corresponding eigenvalues. The 
coefficients c„ are given in terms of the initial wavefunc- 
tion by c n =(4> n (x)\ip(x, 0)). In general the superposition 
(1) may also contain continuum states for which n would 
be a continuous index and the sum would be replaced 
by an integral. Here, however, we concern ourselves with 
only superpositions of bound states assuming negligible 
continuum contributions. Also, one assumes that the ex- 
pansion (I) is strongly weighted around a mean value, 
n. Both the assumptions above are reasonable in, e.g., 
the experimental situation when a localized wavepacket 
is produced using a short laser pulse Q. If one assumes 
that the weighting probabilities |c„| 2 are strongly cen- 
tered around a mean value n, one can expand the energy 
in a Taylor series in n around n as: 



E n = E n + E' n (n -n) + -E%(n - nf 



(2) 



where the primes denote derivatives with respect to n. 
From (2) one can identify the time scales, 



T r j. = 
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(3) 



V(x) = 0, 



, L 

N<2 



(6) 



V(x) = V , \x\>-. 



The system thus has three distinct regions and the so- 
lutions to the Schrodinger equation gives us the energy 
eigenstates and eigenvalues. The even parity solutions 
for the three regions are: 



i(x) — A e^ n cosav, 
(x) — A cos 2a n x, 
(x) — A e' 3 " cos a Tl 



-20 n x 



(7) 



where the eigenvalues are evaluated by solving the tran- 
scendental equation 

a„tana„ = /3„, (8) 

while the odd parity solutions are given by: 



<t>i{x) 



-A e 



0n 



sin a r , 



t >u(x) = A sin2o; n x, 



<I>iii{x) = A e@ n sma„ e 
with eigenvalues given by the transcendental equation 



-2f} n x 



(9) 
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a n cota„ = -p n . (10) 

Here x — x/L where L is the length ol the square well, 

a n = ^JmEL 2 /2Ti 2 and (3 n = ^m{V a - E)L 2 /2h 2 , with 
E the energy (E < Vq) and m the mass of the particle 
and the normalization constant A is given as: 



A = 




(11) 



It is clear from (7) and (9) that the wavefunction has 
nonvanishing values in both the 'classically forbidden' 
regions / and III. The quantum mechanical probabil- 
ity for the particle to be somewhere in regions / and 77/ 
is, therefore, nonvanishing. It is clear, though, that as 
one goes away from the boundaries, the probability den- 
sity decreases rapidly to zero. Thus the particle with 
E < Vq cannot really escape to infinitely long distances 
but stays 'bound' to the well. One can define 



mV L 2 



2hr 



(12) 



as the 'well-strength' |14| . A finite well of well strength e 
would contain a finite number of bound states, N, where 
N ~ — + 1. The time evolution of any given initial 
states can be expressed in terms of these eigenstates and 
eigenvalues, 



N 



(x) exp (- 



%ia 2 T , 



(13) 



where r = t/T rv is the time scaled in terms of the revival 
time T rv ~ AmL /irh of the infinite well potential. The 
expansion (13) contains both odd and even parity states. 
For our purpose we examine the quantum dynamics of an 
initial Gaussian wavepacket with mean position Xq and 
zero mean momentum: 



tp(x, 0) = Bcxp ^ — 



(x - xof 
2a 2 



(14) 



where B is the normalization constant. The time evo- 
lution involves the energy eigenvalues, a n s and the cor- 
responding eigenfunctions which are obtained by solving 
the transcendental equations (8) and (10) numerically. 
A Gaussian wavepacket of zero mean momentum can be 
faithfully constructed by a superposition of these bound 
states ( the sum of the coefficients, |c n | 2 « 1 ). 

Our numerical simulations for the finite well show the 
presence of revivals in the dynamics of the initial Gaus- 
sian wavepacket (14) similar to that seen in the case of 
the infinite well potential for short times. The revival 
times, however, are in general longer than that of the 
infinite well and depend on the well strength (depth), e, 
of the finite well. Barker et al. [[lJJ have shown via a 
first order Taylor series expansion of the transcendental 
equations (8) and (10) for the eigenvalues that the energy 



levels of a finite well of length L and well strength e can 
be approximated as: 



EL = 



TT 2 h 2 n 2 
2mL 2 



(15) 



One can see that this is equivalent to the energy levels 
for an infinite well but with the larger length, 



V = L{1 + - 



(16) 



Corresponding to this approximate expression we can 
thus find an expression for the approximate revival time, 
T' rv , for the finite well in terms of T rv , the revival time 
for the infinite well: 



T' — T 



(17) 



The formula of Barker et al. is more accurate for tightly 
bound eigenstates (deeper wells) than for weakly bound 
states (shallow wells) |14j 

The revival features for the dynamics of any ini- 
tial state can be understood by examining the absolute 
square of the autocorrelation function: 

\A(t)\ 2 = mx,0)i>(x,T))\ 2 =|^|c„| 2 e - i ^|2 (18) 



E 



c B | 2 e- 8fa » T /' r | 2 . 



Fig. 1 shows the square of the autocorrelation function, 
|^4(t)| 2 , for an initial Gaussian wavepacket for three dif- 
ferent 'well-strengths', e, contrasted with that for the in- 
finite potential. One can see that for larger values of 
e, the revival time approaches that of the infinite well 
(T,'^ — > T rv ). The agreement of the actual revival times 
with the formula of Barker et al. gets more accurate as e 
increases, as expected and as shown in Table I. Thus for 
short times the dynamics of the wavepacket in the finite 
well is similar to that in the infinite well with modified 
revival times which depend on the depth of the well. For 
our simulations we work with the values of the parame- 
ters for the initial Gaussian state and the well-strength, 
e such that ^ \c n \ 2 w 1 always. In Fig.l, note that the 
detailed behaviour of the autocorrelation function for the 
infinite well (solid line) shows a symmetry. In contrast, 
the detailed behaviour of the autocorrelation function for 
the finite well (dashed line) shows an asymmetry espe- 
cially around the revival time as is quite evident from 
Fig. 1(a) which corresponds to the shallowest of the three 
wells, (e = 12). As the well gets deeper, the asymmetry 
decreases and as expected, the detailed behaviour begins 
to resemble closer to the infinite well as is obvious from 
the 'deeper' well shown in Fig. 1(c). 
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TABLE I. Comparison of the revival times, T rt)a 
with the revival times estimated from the ap- 
proximate formula of Barker et al, T rv t, for var- 
ious well-strengths, e. All times are scaled by 
T rv = 4mL 2 /irh, the revival time for the infinite square 
well. 



FIG. 1. Square of the autocorrelation func- 
tion, |A(r)| 2 vs the scaled time r for an initial 
Gaussian wavepacket of width L/10 centered ini- 
tially at L/5 for the infinite well (solid-line) and 
finite well (dashed line) of well-strengths: (a) 
e = 12(JV - 8), (b) e = 30 (JV ~ 20), (c)e = 100(JV ~ 64). 
(r = t/T rv ,T rv = AmL 2 /-Kh) 



e 

(well-strength) 


T-rvo, 

(4mL 2 /nh) 


(4mL 2 /nh) 


percentage error 


12 (N ~ 8) 


1.185 


1.174 


0.9% 


30 (N ~ 20) 


1.06777 


1.06876 


.09% 


100 (N ~ 64) 


1.02009 


1.02005 


.0039 % 
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III. SUPERREVIVALS 

In the preceding section we examined the wavepacket 
dynamics and the structure of fractional revivals at short 
times, i.e., times of the order of or close to T', v . Now 
we turn to the behaviour at longer times. A look at 
the square of the autocorrelation function, |v4(t)| 2 , for 
longer times shows that after a few revival cycles, the 
wavepacket ceases to 'reform' faithful to its original form 
at t = 0. The peaks in |^4(r)| 2 signifying full revivals 
go down and then pick up again after a few cycles and 
the wavepacket is once again closer to its original form. 
This new sequence of revivals is characterized by a longer 
revival time (T sr ). This behavior depends, as expected, 
on the depth of the finite well and is more obvious when 
we are dealing with shallow wells, i.e., for smaller values 
of e. For larger depths, e, the usual revivals continue 
for more number of cycles before the peaks in |^4(t)| 2 go 
down and pick up again. The superrevival times, thus, 
get longer and longer with increase in the well depth. 
This is in agreement with our expectation that the be- 
havior in a deeper finite well should approach that for 
the infinite well. It is interesting to note that while the 
usual revival times at short times decrease with increas- 
ing well-strength, e, and approach that of the infinite 
well, T rv = AmL 2 /irh, the superrevival times increase 
with increasing e approaching T sr — oo for the infinite 
well. Since the finite well eigenenergies and eigenfunc- 
tions have no simple analytical form, it is not straight- 
forward to esthave no simple analytical form, it is not 
straightforward to estimate these superrevival times for 
various well-strengths e. Figs. 2, 3 and 4 show the square 
of the autocorrelation function, |^4(r)| 2 , vs r for an initial 
Gaussian wavepacket for two different vvs r for an initial 
Gaussian wavepacket for two different values of the well- 
strength, e. The structure of the plot of the square of the 
autocorrelation function vs the time r is a reflection of 
the initial state and depends on the parameters like the 
width and the mean position of the inon the parameters 
like the width and the mean position of the initial Gaus- 
sian wavepacket. The arrow in these figures indicates the 
first superrevival. We have taken the case of an initial 
wavepacket which has its mean position at the center of 
the well (xq — 0.0) for Figs. 2 andposition at the center 
of the well (xq = 0.0) for Figs. 2 and 3 and (xq = 0.2L) 
for Fig. 4. 
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FIG. 2. Square of the autocorrelation function, 
|j4(t)| 2 vs the scaled time r for an initial Gaussian 
wavepacket of width L/10 and with mean position at 
x = 0.0 for e = 12(N ~ 8); (c) and (d) highlight the de- 
tailed behaviour at shorter times and (e) highlights 
the detailed behaviour near the superrevival time(e) 
highlights the detailed behaviour near the superre- 
vival time(e) highlights t he detailed behaviour near 
the superrevival time. The arrow marks the first su- 
perrevival. The dashed curves correspond to the cor- 
responding plot for the infinite square well. 

One can see that there are many revivals much before 
the characteristic time T' rv (e.g., there are eight revivals 
in Fig. 2(c) and Fig. 3(c) including the one at T' r ^). Aron- 
stein and Stroud JTo| have pointed out this feature in 
their study of wavepacket dynamics in an infinite square 
well potential. They show that for even parity states, 
such as the one where the wavepacket is initially posi- 
tioned at the center of the well (xo = 0.0), there are full 
revivals at multiples of t — T rv /8. Similarly, odd parity 
states show full revivals at multiples of times t — T rv /A. 
These features are obvious in the case of the dynamics in 
the infinite square well potential when one looks at the 
time evolution of an initial state in terms of the even and 
odd parity eigenstates: 
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n(even) 

+ Y \/2r 



e 27rtn r c n sinnTC 



(19) 



c„ cos nnx. 



Here r = t/T rv is the scaled time and x = x/L. For an 
initial even parity eigenstate (ip(x) — ip(—x)), only the 
second summation with the cosines contributes. It can 



be checked that at t = T r , 



this sum can be written as: 



le 27rm ~ T c n cosnirx 



(20) 



i/>(x,t)= Y ^ e 

n(odd) 

= Y^ cos (2mTTx)e-" m ^ m+1 h 

m 

= e~^ip(x,0); m = 0,1,2,.... 



At t — T rv /8 the state is the same as what it was at t = 
except for a constant phase factor. Thus there will be 
full revivals at multiples of t = T rv /8. Similarly, for the 
odd parity states (ip(x) = —ip(—x)) the contribution will 
only be from the first summation containing the sines. 
At t — T rv /i, the time evolved state can then be written 
as: 



ip(x,r)= Y V2e- 2 * m ~ T c n smmrx (21) 



-2-irin 2 T 

v 

n(even) 

^Y, \/2sin (2mwx)e~ 



ijj(x,0); m = 0,1,2, 





Thus full revivals occur here at multiples of t = T rv /A. 
For the case where the initial Gaussian has mean posi- 
tion #0 — 0.2L (a state which does not have a definite 
parity), there is a complete revival only at T rv though 
partial revivals can be seen at multiples of t = T rv /3 
(Fig. 4(c)). These symmetry features have an interest- 
ing analog in the regeneration characteristics of a field 
E(x) of wavelength A p ropagating through a multimode 
planar waveguide of width b [^5|. If the field is sym- 
metric in the transverse dimension (E(x) = E(—x)), its 
regeneration length is L = b 2 /A. An antisymmetric field 
regenerates at a distance 2L and arbitrarantisymmetric 
field regenerates at a distance 2L and arbitrary fields re- 
produce after 8L. In the case of wavepacket dynamics in 
square well potentials, these symmetry aspects are more 
clearly understood for the infinite square well potential. 
However, it is not surprising that we see similar features 
for wavepacket dynamics in the finite well analyzed here. 



FIG. 3. Square of the autocorrelation function, 
|A(t)| 2 vs the scaled time r for an initial Gaussian 
wavepacket of width L/10 and with mean position 
at x = 0.0 for e = 15 (JV ~ 10). (c) and (d) high- 
light the detailed behaviour at shorter times and (e) 
highlights the detailed behaviour near the superre- 
vival time(e) highlights the detailed behaviour near 
the superrevival time(e) highlights t he detailed be- 
haviour near the superrevival time The arrow marks 
the first superrevival. The dashed curves corres- 
The arrow marks the first superrevival. The dashed 
curves correspond to the corresponding plot for the 
infinite square well. 

Also, notice that for the initial Gaussian wavepacket 
centered at xq — 0.0, the autocorrelation function, 
|^4(t)| 2 is never zero (Figs. 2 and 3). On the other hand 
if we examine the case where xo = 0.2L (Fig. 4), there 
will be times when |A(r)| 2 is zero and these times nor- 
mally correspond to those instants when the wavepacket 
appears almost as a mirror reflection of itself and hence 
completely uncorrelated with what it was at t — 0. Lo- 
cal extrema of the autocorrelation function, on the other 
hand, correspond to 'fractional revivals', i.e, when the 
wavepacket appears in a spatially separated superposi- 
tion of replicas of itself (cat-like states). As mentioned 
before, the absence of a simple analytical form for the 
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energy eigenvalues for the finite well makes it difficult to 
find estimates for the superrevival times. The nature of 
these superrevivals and how often they occur would also 
depend on various factors like the well-strength and the 
initial position of the wavepacket. 



time. Although the system described by (19) is not quite 
the same as the finite well potential, one sees certain 
similarities in the autocorrelation function in both cases. 



Comparison with an anharmonic oscillator 

As mentioned above, the lack of an analytical form for 
the energy makes it difficult to estimate the superrevival 
times for the finite well. However, to get some insight, 
the revival and superrevival structure in wavepacket dy- 
namics of the finite well may be compared with the 
wavepacket dynamics in an anharmonic oscillator poten- 
tial where one has the exact analytical expression for the 
energy. We examine the hamiltonian with nonvanishing 
third order nonlinearity: 



H = / Ui(a t a) 2 +/i 2 (a t a) S 



(22) 



Here a and a) correspond to the annihilation and creation 
operators, respectively. This system has been studied by 
Gantsog et al Jll| , though not in the context of quantum 
revivals. To make a comparison with the dynamics of 
the Gaussian wavepacket in the case of the finite well, 
one can study the dynamics of initial coherent states and 
squeezed states for the system described by (19). We 
assume that «2 << Ml- ^ can be easily checked that the 
revival times and superrevival times for the system (19) 
are given by: 



rv = (l + 3n/?)' 

t sr ■ t rv / 



(23) 



l^4(r)| 2 = I 



\ - a 



■ exp 



(-27riTn 2 -27ri/W)| 2 



(24) 



l^)l 2 = l£ 



2yG 



exp 



(s + l)2"n! \s 
2sa 2 



1 



2'KiTn 2 




where t rv = , is the revival time without the third or- 
der term, i.e., with /12 = and = n is the average 
photon number corresponding to the chosen initial state. 
Initial coherent and squeezed states can be expressed in 
terms of the Fock states, \n) , and one can study their 
time evolution. The autocorrelation for an initial coher- 
ent state with amplitiude a can be easily shown to be 



while for an initial squeezed state, the autocorrelation 
function is given by 



Here H n denotes the nth order Hermite polynomial, s is 
the squeeze parameter jlTj and r = t/t rv is the scaled 



FIG. 4. Square of the autocorrelation function, 
|A(t)| 2 vs the scaled time r for an initial Gaussian 
wavepacket of width L/10 and with mean position at 
x = 0.2L for e = 12(N ~ 8). (c) and (d) highlight the 
detailed behaviour at shorter times and (e) highlights 
the detailed behaviour near the superrevival time(e) 
highlights the detailed behaviour near the superre- 
vival time The arrow marks the first superrevival. 
The dashed curves correspond to the corresponding 
plot for the infinite square well. 

The time evolution of an initial squeezed vacuum state 
(a = 0) in the anharmonic oscillator potential (19) can be 
compared with that of the narrow Gaussian wavepacket 
centered at xq = 0.0 in the finite quantum well. Fig. 5 
shows the autocorrelation function for the dynamics of 
an initial squeezed vacuum in the the anharmonic poten- 
tial (19) with squeeze parameter s — 10, and f3 = .002 
for short times as well as for for times close to the su- 
perrevival times. In both these cases one can see that 
the autocorrelation functions never go to zero unlike the 
case when xq 7^ 0.0 for the packet in the finite poten- 
tial well and a ^ for the coherent state/squeezed state 



7 



in the anharmonic potential (results not shown for the 
anharmonic oscillator model). For this 'symmetric case', 
there are many more revivals at short times and super- 
revivals at longer times in both cases in contrast with 
the case when xq ^ 0.0 and a ^ when there are fewer 
revivals and superrevivals. For the anharmonic oscilla- 
tor the autocorrelation function also indicates the pres- 
ence of fractional superrevivals as seen in the finite well 
potential. In particular one can compare the signature 
of the occurence of superrevivals in both potentials by 
comparing the detailed structure of the autocorrelation 
functions at times close to the superrevivals times. In 
both cases one can see that at times nearing the super- 
revival times, the autocorrelation function begins to look 
more and more like the behaviour near r = and near 
the revivals times. Moreover, at exactly these superre- 
vival times the form of the autocorrelation function for 
the finite well appears very similar to that of the infinite 
well and in the corresponding case for the oscillator, the 
autocorrelation function at these times resembles the be- 
haviour for the case when /x 2 = 0. At all other times 
away from the revival and superrevival times, the auto- 
correlation functions look very different. This is clear 
from Figs. 2, 3 and 4 where the detailed behaviour of 
the autocorrelation function for the well at short times 
and near the superrevivals times are highlighted. Thus 
the gross features regarding the nature of revivals and su- 
perrevivals of wavepackets in the finite potential well are 
quite comparable with the dynamics of coherenpotential 
well are quite comparable with the dynamics of coherent 
and squeezed states in a generic anharmonic oscillator 
potential (19). 

IV. CONCLUSION 

In conclusion, we have studied the wavepacket dynam- 
ics in a finite square well potential in the context of quan- 
tum revivals. We have shown that for short times the 
revival patterns are similar to that seen in the infinite 
well potentials but with a modified revival time which is 
now dependent on the depth of the finite well. For deep 
enough wells, the time scales approach closer to that of 
the infinite well. For longer times, the difference in the 
wavepacket dynamics as compared to the infinite well po- 
tential is manifested by the appearance of superrevivals. 
Superrevivals have till now been predicted for the long 
time dynamics of Rydberg wavepackets S . 



(a) m 




FIG. 5. Square of the autocorrelation function, 
|A(r)| 2 vs the scaled time r, (r = t/t rv ,t rv — 2-k%/^,\) 
for an initial squeezed vacuum with squeeze param- 
eter s — 10 and /3 = .002 for the anharmonic oscil- 
lator potential (19). (c) and (d) highlight the de- 
tailed behaviour at shorter times and (e) highlights 
the detailed behaviour near the superrevival time, 
T 3r — T rv /(3 — 500T r „. The arrow marks the first su- 
perrevival. The dashed curves correspond to the cor- 
responding plot for the case when /12 = 

The analysis in this paper predicts that wavepacket 
dynamics in quantum well systems which can be more re- 
alistically modelled as finite square well potentials rather 
than infinite ones, will show superrevivals in addition to 
the usual revivals. 



One of us (GSA) is grateful to W. Schleich for useful 
discussions on the dynamics of wavepackets. 
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